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We consider a finite Fermi-system where the residual interactions create 
a soft mode of the excitation spectrum. Because of the large vibrational 
amplitude, the standard random phase approximation does not work in this 
situation. We develop a regular method for constructing the anharmonic 
potential and illustrate the application of the formalism by a simple model. 


I. INTRODUCTION 


The appearance of a soft collective mode of vibrational nature is quite common in nu¬ 
clei and in other mesoscopic quantum systems, such as atomic clusters. The conventional 
theoretical way to handle this situation |l]-fjj goes along the standard line: mean held (MF) 
— residual interactions — random phase approximation (RPA). The MF determines the 
symmetry of the system around the ground state and corresponding elementary excita¬ 
tions, fermionic quasiparticles. The residual interactions include both coherent effects and 
collision-like processes responsible for the chaotization of motion and lead from a Fermi-gas 
to Fermi-liquid. The formation of coherent modes is described within the framework of the 
RPA where the quanta of those modes are treated as independent quasibosonic excitations 
[2j.|4|] . In the low-lying states, kinematic corrections due to the fact that the quanta are built 
of fermions, as well as high order dynamic effects, can be accounted for perturbatively p,||. 

The RPA-type theories become insufficient when, in some region of the parameters, the 
vibrational frequency oj approaches zero. The vibrational amplitude then grows oc 1 /a/u; 
revealing the instability of the MF. In a finite system, this is not necessarily a vestige of 
a phase transition or sharp restructuring of the system. Rather it might be a signature of 
the failure of the theoretical consideration based on the picture of harmonic vibrations. For 
instance, a low energy of the first excited quadrupole state does not mean that the nucleus 
becomes deformed. In the regime of large amplitude collective motion, we need to reject the 
harmonic approximation and fold a way of calculating the anharmonic effects which cannot 
be here treated as small corrections. 

The estimates [j6|, as well as more detailed unrestricted MF calculations, show that 
in the quadrupole case the effective potential is close to the 7-unstable 0 quartic one, 
~ /3 4 , and the spherical symmetry of the MF still holds but only on average. The popular 
interacting boson model || with the phonon number fixed by a number of fermion pairs 
cannot describe the soft vibrational bands which stretch to very high spins. In this sense 
the phenomenological models |]9j-|T2fl based on a specific form of the collective quadrupole 
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Hamiltonian are more successful. The onset of deformation occurs beyond the point of the 
RPA instability when some previously unoccupied configurations sharply lower their energy 
as a function of deformation and thereby select the equilibrium static value of 7 , usually 
7 = 0. The macroscopic analog of this scenario would be the first order phase transition. 

Below we consider a typical collective soft mode and show a way of constructing the 
effective nonperturbative anharmonic potential for large-amplitude collective motion starting 
with the full microscopic many-body Hamiltonian. We present only the skeleton of the 
formalism and apply it to the Lipkin model known as a testing ground of various theoretical 
approximations (the first version of this approach was published in ||13||). 


II. GENERALIZED DENSITY MATRIX 

It is convenient to use the operator language with the generalized density matrix (GDM) 
as the main tool (Era- We consider a truncated single-particle space of orbitals |1) of full 
dimension fh The GDM is the set of the operators 

R\ 2 = <A a i j (1) 

where a and a' are fermionic operators in the second quantization (in a similar way, one 
can consider Bose-systems), the subscripts 1,2 form a matrix in single-particle space D, and 
each element R\ 2 is an operator in the many-body Hilbert space. These operators generate 
a closed SU(Q) Lie algebra given by the commutation relations 

[L?l 2 j -R3 4] = 4 i? 3 2 — $2 3 R \ 4 ; (2) 

the trace of the GDM in single-particle indices (tr) gives a number operator, tr(i?) = N. 
The GDM is Hermitian in the combined space of single-particle and many-body variables, 
-R12 = -^21- 

The dynamics of the system are governed by a standard Hamiltonian which contains 
one-body and two-body terms, 

H = ^2 e 12«1«2+T X] kl 2; 34 ofll2®3®4 j (3) 

1,2 4 1 . 2 , 3,4 

where e 12 — 621 j Li 2 ; 3 4 = V^* 3 . 2 i, and we assume the antisymmetrized form of the two- 
body interaction. We define a self-consistent held W (similarly to R , an operator in the 
combined space) as a linear functional of the GDM, 

Wia{R) = — X/ ^12; 34 -R 32 = - ^12; 34^4 °3 ■ (4) 

^ 2,3 ^ 2,3 

The hamiltonian in eq. (|3|) can also be written in terms of the GDM. 

The equations of motion for the creation and annihilation fermionic operators are 

[ai, H] = (ei 2 + W 12 ) a 2 , [a{, H} = -'£a\ (e 21 + W 21 ) , (5) 

2 2 

whereas the total GDM ([[]) satishes the nonlinear operator equation 
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[R,H] = [e + W{R},R]. (6) 

These equations are still exact. Here the commutators are understood to act in the combined 
space, for example, 

[w, R] 12 = ^(^13^32 - Ri 3 W 32 y, (7) 
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all elements are many-body operators. 


III. MAPPING ONTO COLLECTIVE SPACE 


Now we make two crucial assumptions: (i) there exists a “collective band” as a set of 
stationary states which are coupled by strong intraband transition amplitudes while the 
transitions to the states of a different nature are weak and can be ignored, or be taken into 
account perturbatively later; (ii) the nomenclature (quantum numbers) of the band states 
can be built with the aid of the operators of collective coordinates a and conjugate momenta 
7 T. These assumptions are fulfilled accurately |I8| for low-lying quadrupole vibrations in 
medium and heavy spherical nuclei where it is known that the quadrupole transitions from 
the ground state are nearly saturated by the first excited 2 + state, which in turn gives 
rise to transitions to the “two-phonon” triplet of states 0 + , 2 + and 4 + , and so on. This 
means that there is a good correspondence between the ideal quadrupole phonon space and 
realistic spectra in spite of the fact that the predictions of the naive model of harmonic 
quadrupole vibrations are badly violated. If so, the observed states can be generated by 
the quadrupole coordinate and momentum operators a 2ll and tt 2 ^ although the collective 
Hamiltonian H(a,ir) might be very far from the harmonic one. 

According to our assumptions, the collective subspace is spanned by the operators a and 
7 r with normal commutation relations 


[a, 7r] = i , 


( 8 ) 


and their high order products. For simplicity we take here scalar quantities; the rotational 
tensor character can be introduced in a straightforward way. The general form of the effective 
Hermitian collective Hamiltonian acting within this subspace is 


m, n 


A( mn ) 
2 mn 


r _ ra _7il 

[a , 7T J + . 


(9) 


Our goal is to derive the unknown c-number coefficients A( mn ) from the microscopic Hamil¬ 
tonian H, eq. (|3|). This can be done by the corresponding mapping of the exact operator 
equations of motion (g). 

We are interested in the matrix elements of the equations of motion between collective 
states. Since the dynamics are saturated in the collective space, we leave as the intermediate 
states in those equations only the states within the band. Then operators R and W can be 
effectively represented by the functions of a and i r similarly to (0), 


^ = E 

m, n 


r (mn) 

2 mn 


\a 


7T 


1 + 5 


w = £ 

m, n 


w (mn) 

2 mn 




( 10 ) 
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The question of mapping is now formulated as a problem of finding a set of numbers hS mn > 
and quantities (matrices in single-particle space) which express the contribu¬ 

tions of specific elementary excitations to a given collective operator. The corresponding 
parts of r( mr d and w^ mn ' > are interrelated by the self-consistency conditions (£|). 

Physical arguments of time-reversal (T) invariance and the possibility of canonical trans¬ 
formations in the collective space, such as shifts and rescalings of collective variables, allow us 
to consider only with even n and start the sum with the harmonic terms, ( mn ) = (02) 

and (20), so that eq. (I) becomes (A^ 02) = 1/B, A^ 20) = C give the mass and force param¬ 
eters of the harmonic part) 


1 2 C 2 A( 3 °) - A^ 12 ) 21 A( 4 °) 4 A (04) 4 A( 22 ) 

H -—n + —a H-— a H-^[a,7r] + H- —a H- —n H-— 


a 2 , 7t 2 1 


+ ... . 
( 11 ) 


Under our assumptions, the full operator equations of motion (|3j) should be satisfied 
inside the band. Therefore we require that in this space 


[K, H + e + W{77}] = 0 . 


( 12 ) 


Commutators involving 7i and e in the above expression are simple since 7i does not con¬ 
tain single-particle variables, whereas e is a c-number matrix in the Hilbert space. The 
commutator of W with the GDM is more complex as both operators act in the combined 
space. 

Below we show the lowest order equations. As seen from ([L2D, it is convenient to introduce 
a self-consistent MF Hamiltonian as 


h = e + W{p} , 


(13) 


and a self-consistent RPA operator L defined 
particle matrix r, 


I5|,[T6| by its action on an arbitrary single- 


Lr = [h,r] + [W{r}, p] . (14) 

In eqs. ( |T5| ) and ( |H| ) we used the ground state single-particle density matrix p = r^ 00) . The 
lowest static part, (mn) = (00), produces a set of the MF equations 

0=[h, P }+i8(°°\ (15) 


where is a correction to the usual Hartree-Fock approximation from higher orders which 
changes the average MF single-particle occupancies (eigenvalues of p) due to the fluctuation 
effects coming from the soft mode |T7[. The next set of equations corresponds to the parts 
linear in a and n operators (T-even and T-odd terms, respectively), 


-iM 2 °) r W=Lr^+i5( 10 \ 


(16) 


, A ( 02 ) r ( 10 ) = Lr( 01 )+^ 01 ). (17) 

These terms are analogous to the RPA although it is not assumed that the occupation 
numbers are 0 and 1. The following three equations in quadratic order are 
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-iA (20 V (11) - iA( 3t V 01 ) 


(18) 


= iLr^ + [w^ 0 \r^]+tS^°\ 


iA (20) r (02) +iA (02 )r ,(20) +iA (12) r (10) = £,.(11) + [^(10) ^(01)] + [^(01) ? r (10)] + ^(11) , 

(19) 


i A (02)r(n) - *A (1 V 01 ) = Hr^ + [w^ 01 \ r^ 01 )] + z^ 02 ). 


( 20 ) 


We limit ourselves here to the fourth order of anharmonicities, i.e. cubic operators in 
equations of motion. The four corresponding equations are 

_ l A (20) r (21) _ zA (30 )r .(ll) _ zA (40) r (01) = I Ar (30) + I ^(2 0)^10)] + I ^(10), ,.(2 0)] + ^(30) ? 


-A\ (20) r (12) 4- *A (02) r (30) - fA (30) r (02) + zA (12) r (20) 4- -A( 22 ) r ( 10 ) = 

2 

Hr (21) + \ [«/ 20 ) ,r (01) ] + [ W ( 11 ),r (1 °)] + I [ w ^\ r ( 20 )] + [w^ 10 \rH + ^ , 

_ iA (20)_(°3) , - A (02) r (21) , 1 A (12) (11) _ * A (22) (01) = 

2 2 

H r ^) + I [ w (° 2 ), r (10) ] + [u/ 1 V (01 )] + [t// 01 >,r (11) ] + 1 [w< 10 >, r (° 2 )] + z^ 12 ) , 

l A (02) r (12) _ i A (12) r (02) + . A (04 )r .(10) = i ^(03) + i [^(0 2), ,.(01)] + I [^(0 1) ? ,.(0 2)] + -^(03) (21) 

2 2 3 2 2 

The higher order corrections 5b4) arise from the commutators [JZ, W] and [JZ, Ti\ , 


s m = I 
2 


n;( 10 ), r< 01 >' 


w (oi) r (io) 


+ ... 


fI A (20),.(13) _ I A (02) (31) , 2.(30) (0 3) _ 1 A (12) (21) , 

V2 2 3 4 


( 22 ) 


Each next term denoted by dots is four orders higher than the previous one. Furthermore, 
the lowest correction due to [JZ, W] is always two orders higher, while the terms from 
[JZ, TL\ are four orders higher, than the similar terms in the left hand side of (|2l|). Their 
contributions become less important || because there the small statistical weight oc 1 /a /O, of 
the collective mode is not sufficiently compensated by the inverse powers of the low frequency 
u. We can note parenthetically that this compensation can occur only in finite systems so 
that the whole approach is tailored for mesoscopic physics. Since the typical estimates for 
the realistic soft modes show the dominance of the quartic anharmonicity, we keep the main 
corrections to the RPA terms 


6 (1 °) 



4 20 >, r< 01 > 

_ 

w V> l>t r ( 20 )' 

A 

U“V (11, 1 

V 


+ 


+ 
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(23) 


5 (01) = l f 

'/V' 01 )' 


W (°D, r (n)' 

+ 

w d°), r (0 2)- 

2 V 


+ 


+ 

■ 




( 24 ) 
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Let us stress here that the method suggested above differs from numerous attempts at 
boson expansion, see || and references therein. We do not map the wave functions from the 
microscopic space to a bosonic one. We also do not map directly the operators of observables. 
We map the equations of motion explicitly truncating the intermediate states. This method 
is regular, does not violate general principles, and, being in fact variational, allows for the 
further improvements by including other intermediate states. 


IV. LIPKIN-MESHKOV-GLICK MODEL 


Using the commonly accepted procedure of testing the validity of many-body approxima¬ 
tion techniques, we apply the method to the two-level Lipkin-Meshkov-Glick (LMG) model 
|L§. The space contains two single-particle levels of energies ±e/2 with a large degeneracy 
0/2 of each of them. We label the O fermionic states by quantum numbers ( crl ), where 
a = ±1 denotes one of the two single-particle levels and Z = 1,2,... ,0/2 distinguishes the 
degenerate states on each orbital. The many-body Hamiltonian of the system is 

H = aa l, i a °,i + \ V a l,i a l, v v i ■ ( 25 ) 

Z <y,l Z a, l,l' 

The special feature of the problem is that the collective dynamics are expressed in terms of 
the quasimomentum operators J±, J z , 

j+ J — Jx T iJy / ) ®+l, / 1, It Jz 7) 'y / ^^< 7 , l l■ (26) 

l Z <7,1 


The Hamiltonian (^5[) can be expressed as 


H - <■/, + -VUl + J 2 _) = <J, + V(Jl - ./=). 


(27) 


The LGM model is ideally suited to our approximate mapping procedure. The SU(2) 
symmetry of the problem can be also combined with particle-hole symmetry, which allows 
us to limit the consideration to the cases with the particle number N < Vt/ 2, and discrete 
symmetries (we can take V > 0 without loss of generality). For the unperturbed system, 
V = 0, the ground state with all N particles on the lower level belongs to the largest 
representation J = J z = N/2 with J 2 = [(N/ 2) + l](iV/2), and then J + is an operator 
that creates a collective state. In this model the collective degrees of freedom are decoupled 
exactly, and we need to reproduce the equations of motion 


[ Jx■ H\ ieJy iV^Jy^ T 2 ]_)_, \Jyi JJ\ i^Jx d-]. , 


(28) 


[Jz, H} = 2 iV[J x , Jy} + , 


(29) 


in the mapped space of collective variables a , 7r, with a collective anharmonic Hamilto¬ 
nian (|TT|). The kinematic constraints, analogous to eq. (|^, arise from the mapping of the 
quasimomentum algebra onto a Heisenberg algebra of a and tt. These constraints can be 
accounted for by the Holstein-Primakoff transformation 
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J+ = Jl = A^2J -A^A = \j2J + 1 — A^A A\ J z = -J + A'A , (30) 

where operators A and A^ are bosonic annihilation and creation operators, [A, .47] = 1. 

The RPA corresponds to keeping only quadratic terms in A and A^ which leads to the 
RPA Hamiltonian 

H rpa = -e (j + i) + e -{A^A + AA f ) + jVl6J 2 -l ((A^) 2 + . ( 3 l) 

The diagonalization of ([H]) results in the harmonic approximation with an RPA frequency 

2Q 

<4pa = e 2 - V 2 (4 J' 2 - 2 ) • (32) 

The instability point, V 2 ~ e 2 /(4J 2 ), exists in all J-subspaces emerging first for the largest 
J with a greater degree of collectivity. 

The collective coordinate and momentum can be introduced with the aid of the canonical 
transformation 

A 7 = — y= (—iua + vir ), uv = —1. (33) 

The LMG model has only even order anharmonicities, and in our choice of expansion the 
correction to the n-th order will come from the (n + 2)-th order in a and n. An expansion 
up to the sixth order retaining only quadratic and quartic terms is necessary for determining 
the effective quartic Hamiltonian. The appropriate choice of u and v as 


A = (iua + vi r), 

V2 


u= e + V 


16 J 2 + 8J - 1 


1/2 


v = —, 
u 


2V16 J 2 - 1 ) 

sets a scale of the collective Hamiltonian at B = 1 and the parameters in eq. Q as 

r2 (16 J 2 + 8 J — l) 2 


co z = C = e z 


V 2 


4(16 J 2 - 1) 


4 V 2 J 2 , 


(34) 


(35) 


A ( 4 o) 


/ 32J 3 — 2 J — l\ 
V (16J 2 - I ) 3 / 2 J 


u 


Vu\ 


x m 


—2V 


( 32J 3 - 2J- 1\ 

V ( 16 J 2 - I ) 3 / 2 J 


v 4 ps — V v 4 . 

(36) 


At the instability point of u —> 0, assuming that J 1, e 3> R, we obtain an approximate 
collective Hamiltonian (V > 0) 


H 


2 


+ 4e 4 V a 4 . 


(37) 


The negative 7r 4 term in the collective Hamiltonian is very small in the vicinity of the 
instability point in contrast to the quartic potential a 4 which has large matrix elements 
because of the large amplitude of collective motion. This is a typical situation which emerges 
with a suitable choice of collective coordinates (T-even) and collective momenta (T-odd). 
The next order terms and, in general, coupling to non-collective space, will correct the 
behavior of the n 4 term but this is not important for low-lying states. 
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FIG. 1. Parameters of the collective Hamiltonian (harmonic frequency and quartic parameters) 
for the LGM model with N = 20 particles in the largest representation J = 10 as a function of the 
interaction strength V (e = 1), upper panel (note the logarithmic scale with the strong growth of 
A^ 40 ) at the point of the RPA instability); excitation energies of the first five states for the same 
case of the LMG, exact solution, solid lines; RPA solution, dotted lines; anharmonic oscillator 
solution with the 7r 4 term ignored, dashed lines, lower panel. 

In Fig. [I] we show the behavior of the harmonic term (cu) and quartic corrections A^ 40 ' 
and A^ 04 ) in the dimensionless normalization, upper panel, and present a comparison of the 
exact LMG model spectrum (solid lines), RPA solution (dotted lines), and an improved 
anharmonic oscillator solution with the ignored divergent 7r 4 part (at large V it should 
be included along with the high order coordinate terms), lower panel. The anharmonic 
effects produce a dramatic improvement as compared to the RPA. At the point of the RPA 
instability, the contribution from the large quartic potential restores the stability. As the 
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interaction strength V increases, the effective potential 

U(a) = -cu 2 a 2 « 4 (38) 

evolves, Fig. 2, from the harmonic oscillator to a pure quartic oscillator at the instability 
point, and to the “Mexican hat” potential with two minima. In the last limit the lowest 
states of opposite parities located in the minima become degenerate as clearly seen in Fig. 2 
(“chiral symmetry”). Contrary to the macroscopic second order phase transition, the higher 
states are located above the barrier and feel only the main quartic potential. A similar 
phenomenon should exist in soft nuclei beyond the RPA instability point when only the 
lowest states are influenced by the presence of the minima in the /3 coordinate; however, 
there the minima are connected along the 7 coordinate which is absent in the LGM. 



FIG. 2. The spectrum of the lowest levels in the model of Fig. 1, dashed lines, and the shape of 
the effective potential U(a ) for (a) V = 0, the harmonic limit; (b) V = 0.05, pure quartic potential 
at the RPA instability point; (c) V = 0.08, beyond the instability point, the splitting of the levels 
of opposite parity (symmetry in J z in the limit of V J 2 3> e) below the barrier decreases as V 
increases. 


V. CONCLUSION 

We have suggested an alternative approach to the construction of a collective Hamiltonian 
for large amplitude collective motion in a finite Fermi system in the presence of a soft 
vibrational mode. In such a situation, the RPA is insufficient as near the RPA instability 
the anharmonic effects dominate. The advantages compared to conventional techniques, such 
as the generator coordinate method, are related to a fully quantum consideration which does 
not require the derivation of an approximate classical Hamiltonian with the subsequent ill- 
defined procedure of requantization. It differs from the approaches utilizing various versions 
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of the boson expansion in the variational character of the formalism. We can vary the 
collective space assuming the saturation of the exact operator equations of motion within a 
part of total Hilbert space. With our techniques we also avoid the slow convergence problem 
of the naive boson expansion. 

Of course, the illustrative example of the Lipkin model is perfectly suited to our goal 
since the collective modes of this model are fully decoupled. However it emphasizes the pre¬ 
dominance of the quartic anharmonicity near the RPA instability. Because of the convenient 
operator distinction between the coordinate and momentum parts, we concentrate the most 
important anharmonic effects in the quartic potential which has large matrix elements in 
the dangerous region of interest. 

In realistic cases, the collective space is not decoupled completely. The effects of coupling 
to noncollective states lead to the spreading of the collective strength and the chaotization 
of motion in the region of high level density. To treat this situation as well, we can include 
the matrix elements of the GDM connecting the collective band with incoherent states. One 
promising approach would be to consider these states on average, making the random phase 
assumption on a new level of treatment. This would introduce an effective background for 
the collective mode to describe its spreading and damping. 

The authors are thankful to D. Mulhall for assistance and criticism; they acknowledge support 
from the USA National Science Foundation. 
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